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We analyzed the influence of static gravitational field on the vacuum and proposed the concept of inhomoge- 
neous vacuum. According to the observational result of the light deflection in solar gravitational field as well as 
the corresponding Fermat's principle in the general relativity, we derived an analytical expression of the refrac- 
tive index of vacuum in a static gravitational field. We found that the deviation of the vacuum refractive index 
is composed of two parts: one is caused by the time dilation effect, the other is caused by the length contraction 
effect. As an application, we simulated the effect of the gravitational lensing through computer programming 
and found that the missing central imaging could be interpreted in a reasonable way. 
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I. INTRODUCTION 

Vacuum is usually considered as "homogeneous" and 
"isotropic", i.e., vacuum does not differ from place to place, 
and the refractive index of vacuum is always equal to 1. 
However, the recent theoretical and experimental progresses 
demonstrate that such concept of vacuum turns out to be in- 
appropriate when there are matters or fields within finite dis- 
tance. For example, the vacuum inside a microcavity is modi- 
fied due to the existence of the cavity mirrors, which will alter 
the zero-point energy inside the cavity and cause an attrac- 
tive force between the two mirrors known as Casimir effect 
0112], which has been verified experimentally IH0]. A sec- 
ond example is that, under the influence of electromagnetic 
field, vacuum can be polarized, which has led to astonishingly 
precise agreement between predicted and observed values of 
the electron magnetic moment and Lamb shift, and may in- 
fluence the motion of photons 01 . Dupays et al. |@] studied 
the propagation of light in the neighborhood of magnetized 
neutron stars. They pointed out that the light emitted by back- 
ground astronomical objects will be deviated due to the optical 
properties of quantum vacuum in the presence of a magnetic 
field. Also in |7|], Rikken and Rizzo considered the anisotropy 
of the optical properties of the vacuum when a static mag- 
netic field Bo and a static electric field Eo are simultaneously 
applied perpendicular to the direction of light propagation. 
They predicted that magnetoelectric birefringence will occur 
in vacuum under such conditions. They also demonstrated 
that the propagation of light in vacuum becomes anisotropic 
with the anisotropy in the refractive index being proportional 
to B x E . 

The facts that the propagation of light in vacuum can be 
modified by applying electromagnetic fields to the vacuum 
implies that the vacuum is actually a special kind of opti- 
cal medium (51 [60 . This is similar to the Kerr electro-optic 
effect and the Faraday magneto-optic effect in nonlinear di- 
electric medium. This similarity between the vacuum and the 
dielectric medium implies that vacuum must also have its in- 
ner structure, which could be influenced by matter or fields as 



well. Actually, the structure of quantum vacuum has already 
been investigated in quite a number of papers MM> FToh. 

In this paper, with the analysis of the influence of static 
gravitational field on the vacuum, we put forward a new con- 
cept that the curved spacetime around a certain matter can be 
treated as an optical medium with a graded refractive index. 
We suggest that the so-called curved spacetime is a reflec- 
tion of the vacuum inhomogeneity caused by the influence of 
gravitational matter. Based on this idea, the refractive index 
of vacuum is derived. We will also apply this concept to un- 
puzzle the problem of the central image missing in almost all 
the observed cases of gravitational lensing 1111m . 



II. THE DEVIATION OF THE VACUUM REFRACTIVE 
INDEX 



According to the astronomical observation, the light prop- 
agating through a space with a celestial body nearby will be 
deflected. It can be interpreted with the curved spacetime in 
general relativity. As a matter of fact, it can also be inter- 
preted with the assumption that the vacuum around matter is 
inhomogeneous with refractive index deviated from 1 . Here 
we put forward a theoretical model to describe the refractive 
index profile based on the Fermat's principle for the propaga- 
tion of light in a static gravitational field, which was given by 
Landau and Lif shitz 11 1 211 : 
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where dl is the local length element passed by light and mea- 
sured by the observer at position r in the gravitational field, r 
is the distance from this element of light to the center of grav- 
itational matter M, goo is a component of the metric tensor 



1/2. 
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gfiv, #oo ' ^1 corresponds to an element of optical path length, 
goo" 1/2 = dt/dr, where dr represents the time interval mea- 
sured by the local observer for a light ray passing through the 
length dl, while dt is the corresponding time measured by the 
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observer at infinity. Eq.(l) could then be rewritten as 



goo 



■ l/2 dl 



J dr 

C dt dl 
J dr ds 



nds = 0, 



(2) 



where ds is the length element measured by the observer at 
infinity, corresponding to the local length dl. 

Eq.(2) shows that if we set the scale of length and time at 
infinity as a standard scale for the whole gravitational space 
and time, the propagation of light then satisfies the standard 
representation of Fermat's principle, with the space — actu- 
ally the vacuum — possessing a refractive index given by 



dt dl 
n = — — = n\n 2 . 
dr ds 



(3) 



The factor n\ of the refractive index relating to the time 
transformation effect dt/dr can be derived from the Newto- 
nian attraction, which contributes partially to the deflection 
of light. Considering a photon of relativistic mass at the 
infinity moving down to position r, the work done on to the 
photon by the Newtonian gravity is 



GMm 



-dr = d(mc ), 



(4) 



where G is the gravitational constant, c is the velocity of light, 
M is the mass of a star (say the Sun), r is the distance to the 
center of the star. Integrating Eq.(4) gives 



(5) 



where m r is the relativistic mass of the photon at position r. 

Since the photon energy is E = hv = mc 2 , where h is the 
Planck constant, v is the photon frequency, then we have m = 
hv/c 2 . Substituting it into Eq.(5) gives 



GM 



It is just the frequency shift caused by the gravitational force, 
which reflects that a clock in a gravitational field runs slower 
than that far away from the gravitational center. That is 



dr- 



GM 

" ^ dt, 



(7) 



where dr denotes the time measured by a clock at position 
r, dt is the converted time of dr, i.e., the time measured by 
the clock at infinity. This relation indicates that, if the length 
scale is the same, i.e., dl = ds, when an observer at position r 
reports a light velocity c\ - dl/dr, it should be converted by 
the observer at infinity into 



c, = 



ds 
dt 



dl 



GM 
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FIG. 1 : Light deflection caused by a graded refractive index. 



This change of light velocity will certainly bring a deflection 
to the light propagation. The corresponding refractive index 
is 



dt 



Yl\ - — e rc Z 

c, dr 



(9) 



Let us now consider the deflection angle caused by this 
graded refractive index. In Fig.l, the curve AP represents the 
light ray, ft is the angle between the position vector r and the 
tangent at the point P on the ray, <p is the deflection angle of 
light. Since the refractive index shown in Eq.(9) has spheri- 
cal symmetry, i.e., depends only on the distance r for a given 
mass M, according to the Fermat's principle 



nds = 0, 



(10) 



where ds = dr + (ra) 2 , a = da/dr, n = n(r), we have the 
corresponding Lagrangian function 



L(a, a\ r) = n(r) y/T+ (ra) 2 . 



Using the Lagrangian equation 



(6) we get |U3| 



or 



dr da da 



nr sin J3 = constant, 



nr sin/3 = n^ro, 



(11) 



(12) 



(13) 



(14) 



where r and no represent the radius and refractive index at the 
nearest point A respectively. 
Since 



tan [3 = 



rda 
dr 



associating with Eq.(14) reaches 

dr 



da 
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(15) 
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FIG. 2: Light deflection in solar gravitational field. 



For a light ray passing by the Sun as shown in Fig. 2, the 
total angular displacement of the radius vector r reads 



Aa 



(17) 



where ro represents the nearest distance to the center of the 
Sun. Because the gravitational field of the Sun is a weak field, 
the value of GM/rc 2 is quite small, so substituting Eq.(9) into 
Eq.(17) gives a solution of first order approximation 



Aa 



71 + 



2GM 

r c 2 ' 



(18) 



Then the total deflection angle of light caused by the refrac- 
tive index n\ in solar gravitational field is 



A(f\ = Aa - 7i = 



2GM 

r c 2 ' 



(19) 



In fact, this result was obtained early in 1911 by Einstein 
full , who also investigated the effect of red shift and the cor- 
responding slowing down of the light velocity in gravitational 
field and then figured out the light deflection as shown in 
Eq.(19) with the use of Huygens' principle. Since the ac- 
tual total deflection angle of light propagation calculated by 
the general relativity tlHQii and measured by the astronomi- 
cal observation ifnll is twice that value, we then come to know 
that the length transformation effect dl/ds in Eq.(2) must have 
the same relation as that of the time transformation effect 
dt/dr expressed in Eq.(9), namely 



dl QM 

ds 



(20) 



This relation indicates that a ruler in a gravitational field is 
shorter than that far away from the gravitational center. So 
when an observer at position r reports a length dl, it should be 
converted by the observer at infinity into 



ds 



GM 

e «2 dl. 



(21) 



For a light passing through a length dl, if the time scale is 
the same, i.e., dr - dt , the light velocity = dl/dr reported 
by the observer at position r should then be converted by the 
observer at infinity into 



ds 
dt 



_ GM 

e ~^ dl 



dr 



(22) 
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FIG. 3: The dependence of the vacuum refractive index n on the 
distance r. 



This change of light velocity will also bring a deflection to 
the light propagation. The refractive index corresponding to 
this kind of deflection is 



C2 



dl 



n 2 = — = — = 

c 2 ds 



which also causes a deflection angle of light 



Acp 2 



2GM 

r c 2 ' 



(23) 



(24) 



Therefore, the total deflection angle of light in solar gravita- 
tional field is 



A(f - Aip\ + A(p2 



4GM 

r c 2 ' 



(25) 



The above result shows that, if the two refraction effects are 
considered simultaneously, then the gravitational space — ac- 
tually the vacuum in the gravitational field — can be regarded 
as an optical medium with a total refractive index given by 



2GM 
e rc 2 . 



(26) 



n is composed of two factors: n\ — related with the time trans- 
formation or "curved time"; n 2 — related with the space trans- 
formation or "curved space". So the curved spacetime of gen- 
eral relativity is reflected in the synthesized refractive index 
n, which is also a reflection of the inhomogeneity of the vac- 
uum, showing that the vacuum near the matter is influenced 
more than that far away from the matter. 

The above expression of n shows that the refractive index 
of the vacuum at the infinity from the gravitational matter is 
1, i.e., the usual refractive index of vacuum. The closer of the 
position to the center of matter M, the higher the refractive 
index of the vacuum. The relation between n and r is depicted 
in Fig. 3, where 2GM/c 2 is taken as the unit of r. For exam- 
ples, the corresponding radii for the surface of the Sun in solar 
gravitational field and the surface of the Earth in earth gravita- 
tional field are 2.36 x 10 5 and 7.20 x 10 8 respectively — both 
are far beyond the r-axis illustrated in Fig. 3. 

The deviation of the vacuum refractive index from the usual 
value 1 is given by 



An = n - I = e ^ 2 - 1 . 



(27) 



FIG. 4: The two images I\, I 2 of a gravitational lens. 



FIG. 5: A ray tracing result for the central imaging. 



In weak field it becomes 



distribution 



An : 



2GM 



(28) 



In order to provide the readers with a quantity impression, 
let us give two examples. For the solar gravitational field 
(Af = 1.99 x I0 30 kg), the deviation of n on the surface of 
the Sun (r = 6.96 x 10 8 m) is 4.24 x 10" 6 . For the earth grav- 
itational field (Af = 5.98 x 10 24 &g), the deviation of n on the 
surface of the Earth (r = 6.38 x 10 6 m) is only 1.39 x 10" 9 , 
which is so small that it can hardly be observed in usual ex- 
periments. Nevertheless, for a massive celestial body such as 
a heavy star, a galaxy or a cluster of galaxies, the deviation is 
not only observable, but also important and useful in gravita- 
tional astronomy. 



III. APPLICATIONS 



p=p [i-(-n 



(30) 



where po is the central density of the lens, < r < R, k > 0. 
The density p decreases with the distance r from the center of 
mass; the decreasing varies with the parameter k. This model 
gives the distribution of gravitational potential as 



P = -47ip G 



R 3 



Pi = -4np G{ 



2(2 + k) 



-R l 



[i 



1 



3(3 + k) r 



(-)V},(31) 



6 (2 + k)0 + kyR 



for outside (r > Rq) and inside (r < Ro) the gravitational lens 
respectively. 

The vacuum refractive index profile outside and inside the 
gravitational lens then reads 



The deflection of light by massive bodies leads to the effect 
of gravitational lensing. Formerly, this effect should be calcu- 
lated complicatedly with the general relativity iflUl . Once we 
have introduced the concept of graded vacuum refractive in- 
dex and obtained its relation with mass Af and position r, the 
problem of gravitational lensing could then be treated easily 
with the conventional optical method. 

Considering a source S and a lens L of mass Af, the light 
emitted from S is bent due to the gravitational field of the 
lens. The bent light could be figured out through Eq.(13) and 
Eq.(26). Drawing the extension line of the light from the ob- 
server 0, the apparent (observed) position of the source image 
I could then be found out. The result is shown in Fig. 4. 

This method could also be applied in studying the central 
imaging. In doing this, the vacuum refractive index profile 
inside the lensing body should be considered as well. 

Noticing that Eq.(26) could be virtually rewritten as 



n-e c 2 



(29) 



where P r represents the gravitational potential at position r 
from the center of the lens. 

As a model for discussion, we suppose a lens (for example, 
a galaxy or a cluster of galaxies) of radius R with a density 



exp 
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Fig. 5 shows a ray tracing result for the imaging of a gravita- 
tional lens with the above described vacuum refractive index 
profile. In the figure, only three paths of ray (the three thick 
lines) could pass through the observer 0, forming the upper, 
lower and central images respectively. From the figure, we 
find that, under the same conditions, the larger the distance 
OL from the observer to the lens, the closer the central imag- 
ing light to the center of the lens. If the source S and the 
observer O are counterchanged, it could also be known from 
the figure that, the larger the distance S L from the source to 
the lens, the closer the central imaging light to the center of 
the lens. In addition, through the change of the lens mass 
Af = ^nR 3 pok/(3 + k), we also find that, when the mass Af 
increases, the distance from the central imaging light to the 
center of the lens decreases (Fig. 6, where the mass ratio of the 
lenses corresponding to the four central imaging rays from 
bottom to top is 2 : 3 : 4 : 5 ). 

For the actual condition of gravitational imaging, the dis- 
tances OL, SL and the mass Af are all astronomical figures; 
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IV. CONCLUSIONS 



FIG. 6: Tracing the central imaging rays for lenses of different mass. 

therefore, the light of central imaging is extremely close to 
the center of the lens. However, for a lensing body with a 
density increasing towards the center, it is possible that there 
are barrier matters near the center which will destroy the for- 
mation of the central image. Besides, the relatively longer 
inner path of the central imaging light adds the possibility of 
light being held back by the lens matters on the way. These 
and some other factors such as the relative faintness of the 
central imaging light and the possibly higher brightness of the 
lens core itself, all decrease the possibility of central imaging 
being actually observed. This analysis is firmly supported by 
the fact that the number of observed images is not "odd" as 
expected by the existed theories but "even" in almost all cases 
of gravitational lensing II 1 IN - 



We have proposed the concept of inhomogeneous vacuum 
with graded refractive index based on the analysis of the influ- 
ence of static gravitational field on the vacuum, we derived the 
expression of this refractive index analytically. By using this 
expression, we investigated the effect of gravitational lensing 
in a conventional optical way and provided a reasonable inter- 
pretation for the problem of central image missing. 

The result indicates that, the concept of inhomogeneous 
vacuum is mathematically equivalent to the curved spacetime 
in the general relativity; therefore, an effective and convenient 
alternative method (i.e., optical method) could be established 
to solve the so complicated problems in gravitational astron- 
omy. Physically, under such point of view, the motion of light 
in gravitational space is a motion of light wave in a quantum 
vacuum with graded refractive index. And as we know that, 
in conventional optics, the Fermat's principle says that the op- 
tical path between two given points is an extremum. This is 
also equivalent to the theorem in the general relativity that a 
particle always moves along a geodesic line in a curved space- 
time. 
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